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Abstract. We show that the chiral partition function of two-dimensional Yang- Mills theory 
on the sphere can be mapped to the partition function of the homogeneous six-vertex model 
with domain wall boundary conditions in the ferroelectric phase. A discrete matrix model 
description in both cases is given by the Meixner ensemble, leading to a representation in terms 
of a stochastic growth model. We show that the partition function is a particular case of the 
' z-measure on the set of Young diagrams, yielding a unitary matrix model for chiral Yang- Mills 

theory on and the identification of the partition function as a tau-function of the Painleve V 
' equation. We describe the role played by generalized non-chiral Yang-Mills theory on in 

(N ■ relating the Meixner matrix model to the Toda chain hierarchy encompassing the integrability 

• of the six-vertex model. We also argue that the thermodynamic behaviour of the six-vertex 

model in the disordered and antiferroelectric phases are captured by particular g-deformations 
' of two-dimensional Yang-Mills theory on the sphere. 

cn 



1. Introduction and summary of results 
P-i- 

Qj , 1.1. Gauge theories and statistical mechanics. The relationships between gauge theories 

and exactly solvable models of statistical mechanics has been a subject of intense activity for 
many years, in diverse dimensions and contexts. Such connections have the potential to unveil 

■ the dynamical reasons as to why only particular selected classes of systems are integrable. In 
• recent years this study has also shed light on particular tractable sectors of gauge theories. The 
, most distinctive example is the study of the planar limit of = 4 maximally supersymmetric 
[ Yang-Mills theory in terms of integrable spin chains [H [2] (see [3] for a recent review). In 

^ ' the context of the AdS/CFT correspondence, the integrability structures on both the string 

■ theory and gauge theory sides have been identified and matched. In this paper we present some 
new connections between certain integrable lattice models of statistical mechanics and two- 
dimensional Yang-Mills theory; this theory has a long history as an exactly solvable quantum 
gauge theory with deep connections to one-dimensional integrable systems of Calogero-type, 
string theory, topological field theory, and the geometry of moduli spaces. 



X 



1.2. Two-dimensional Yang-Mills theory. In this paper we study quantum Yang-Mills the- 
ory with gauge group SU{N) on an oriented closed Riemann surface E of genus h and unit area 
form dfj, [1]. The action is 

(1.1) Sym = -^ [ dfiTvF^ , 

where gs plays the role of the coupling constant, F is the field strength of a matrix gauge 
connection, and Tr is the trace in the fundamental representation of SU{N). It was Migdal's 
idea to utilize a lattice regularization of the gauge theory, which relies on a triangulation of the 
two-dimensional manifold S with group matrices situated along the edges [5] . The path integral 
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is then approximated by the finite-dimensional unitary matrix integral 

(1.2) Zm = J Yl dUe II Zp[Up] , 

edges I plaquettes P 

where dUi denotes Haar measure on SU{N) and the holonomy Up = H^eP the ordered 

product of group matrices along the links of a given plaquette. The local factor Zp [U p] is a 
suitable gauge invariant lattice weight that converges in the continuum limit to the Boltzmann 
weight for the Yang- Mills action (jl.ip . 

There are two common choices for the lattice weight Zp [Up], involving the Wilson action and 
the heat kernel action. The latter action has many interesting features [6J and is the usual choice 
in two-dimensional Yang-Mills theory [4J. It leads to the well-known group theory expansion of 
the partition function [5l [7] 

(1.3) Zu = Yl (dimA)2"2/« ( _ ^ 

A 

where the sum runs through all isomorphism classes A of irreducible representations of the 
SU{N) gauge group, dim A is the dimension of the representation A, and C2(A) is the quadratic 
Casimir invariant of A. 

In the lattice approximation one generally expects that in the limit of a very fine triangulation 
we get a theory that converges to the continuum gauge theory. In the present case, this approach 
is in fact much more powerful due to an invariance property of the partition function (jl.3p under 
subdivision of plaquettes of the lattice, which is specific to two dimensions. This property implies 
that the lattice computation is independent of the chosen triangulation of S, and hence is exact 
on an arbitrarily large lattice [5l [8]. This opens up the possibility that there may exist an 
integrable lattice model of two-dimensional statistical mechanics that exactly describes, or is 
equivalent to, two-dimensional Yang-Mills theory in the continuum. In this paper we show that 
this is indeed the case for the chiral sector of the gauge theory, which is defined as follows. 

The sum in ()1.3p runs through all irreducible representations of SU{N), but one can also 
restrict the sum to a subclass of representations. In the large N limit, any representation A of 
SU{N) can be expressed uniquely [9] in terms of coupled representations A = A+0 A_, defined 
to be the largest irreducible representation in the decomposition of the tensor product A+ (8) A_ , 
such that the Young tableau for A is given by joining a chiral tableau A+ to an antichiral tableau 
A_. The number of boxes in the Young tableaux corresponding to X± are understood as being 
small compared to A^. Then the Hilbert space ^ym^'* class functions on SU{N) factorizes 
for large into the coupled tensor product of two sectors [9l [10] 

(1.4) lim n^'i}''^ = n+^n- . 

N^oo 

The chiral Hilbert space ^+ consists of states corresponding to "small" representations A+ of 
SU{oo) in which the number of Young tableau boxes is an arbitrary but finite non-negative 
integer, while the antichiral Hilbert space Ti- consists of states corresponding to conjugates of 
small representations. 

The chiral partition function Zym(S, SU{N)) on an arbitrary surface S is defined by keeping 
only the states of in the large N Hilbert space ()1.4p . This definition also makes sense in 
the g-deformation of the gauge theory (TTJ [T^] that we additionally consider, whose heat kernel 
expansion is given by substituting the dimensions dim A of SU{N) representations with their 
quantum dimensions dim^ A in ()1.3p . In this case the chiral/antichiral factorization is important 
in the context of the OSV conjecture in topological string theory [13j. 
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1.3. Six-vertex model with domain wall boundary conditions. In [8J Witten indicated 
a relationship between two-dimensional Yang-Mills theory and IRF models, in the presence of 
Wilson loops and in a rather generic way. He showed that the lattice gauge theory description 
of Wilson line correlators could be expressed as a lattice statistical mechanics formula similar 
to that of an IRF model. In this paper we explore a different connection with integrable lattice 
models which are dual to IRF models, i.e. six-vertex models. We show that the partition 
function ("^^ ' 'S'C/(A^)) for the chiral sector of Yang-Mills theory on the two-sphere S = 
S*^ can be mapped to that of the six-vertex model with domain wall boundary conditions in 
its ferroelectric regime [T^. The results of this paper give gauge theory derivations of these 
lattice models which elucidate further integrability properties on both the gauge theory and 
statistical mechanics sides. The gauge theories may also provide computationally useful means 
for exploring various aspects and for understanding the origins of these exactly solvable models. 
Previous correspondences between non-chiral two-dimensional Yang-Mills theory and certain 
one-dimensional integrable systems are found in \15\ \W[ [T7] . 

The six- vertex model is a two-dimensional exactly solvable lattice statistical mechanics model, 
introduced by Lieb and Sutherland [El UHl [20] , in which local states are associated with edges 
of an X square lattice and local statistical weights are assigned to its vertices. The states 
each take two values that are usually denoted as arrows along the edge. Since the lattice is 
square, there are in principle 16 possible arrow configurations around each vertex, but most of 
them are chosen to have zero weight in such a way that only six configurations are allowed with 
equal numbers of incoming and outgoing arrows. The partition function of the six-vertex model 
is then 

(1.5) z^= E 

arrow i=l 
configurations cr 

where Wi, i = 1,...,6 are weights associated to each possible vertex state and Ni{a) is the 
number of vertices of type i in the configuration a. 

The six-vertex model suffers from an intricate dependence on boundary conditions, due to the 
constraints imposed by arrow conservation. In particular, the free energy computed with domain 
wall boundary conditions is different from that computed with periodic boundary conditions, 
even in the infinite volume limit |2 H [22 1 [23]. With domain wall boundary conditions, all arrows 
on the left and right boundaries are outgoing, while on the top and bottom boundaries all arrows 
are incoming. In addition to demonstrating the dependence of thermodynamic quantities on the 
boundary conditions [241 114] . there is currently much interest in this model due to its deep 
connections with several problems in algebraic combinatorics. In particular, Kuperberg studied 
its partition function, and used it to give a direct and transparent proof of the alternating sign 
matrix conjecture [251126]. There is also a direct connection with the study of domino tilings [24]; 
in particular, the free fermion line of the six-vertex model with domain wall boundary conditions 
is related to the domino tiling problem for the Aztec diamond |27] . 

Recall the exact self-similarity property of two-dimensional lattice gauge theory [5l [8]. We 
will show that the precise macroscopic two-dimensional lattice model that can be mapped to 
SU{N) chiral Yang-Mills theory on S*^ is the six-vertex model with domain wall boundary 
conditions in the ferroelectric phase, in the limit introduced in [13]. In the six- vertex model the 
rank of the gauge group A^ corresponds to the size of the N x N square lattice; the mapping 
between coupling parameters is given in ()3.8p below. We will also argue that the six-vertex model 
description beyond this approximation constitutes an interesting special case of generalized two- 
dimensional Yang-Mills theory. This gives a statistical mechanics interpretation of the chiral 
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sector of Yang-Mills theory, and also an interpretation of the partition function in terms of 
the normalization constant of a certain stochastic process. Conversely, the six-vertex model in 
this ordered phase and with these boundary conditions can be described as a topological gauge 
theory on a single plaquette. 

1.4. Unitary matrix models and Painleve transcendents. The partition function of two- 
dimensional Yang-Mills theory based on the heat kernel action ()1.3p . i.e. as a sum over ir- 
reducible SU{N) representations A, can be rewritten in the linearized chiral case in terms of 
the normalization constant of the Schur measure [28]. This is especially relevant for the chiral 
sector, because in this case it is straightforward to write the partition function as a Toeplitz 
or Fredholm determinant [28]. This determinant expression leads to connections with unitary 
matrix models and with integrable hierarchies. 

The first physical model whose correlation functions were expressed as a Toeplitz determi- 
nant was the two-dimensional Ising model, arguably the simplest of all integrable systems. The 
seminal work [291 [30] established that the two-point correlation functions of spin and disorder 
fields can be expressed in terms of a solution to the Painleve III equation. This result was 
extended and formalized by the Kyoto school in [3T]. In particular, they connected the the- 
ory of isomonodromy preserving deformations of linear differential equations with the n-point 
correlation functions of the two-dimensional Ising model, and also related the reduced density 
matrix of the impenetrable Bose gas model with the Painleve V transcendent. Other models, 
possessing a free fermion region, also have correlation functions that solve non-linear differential 
equations [32l[33]. In this approach, a Fredholm determinant representation of the correlators 
is crucial. 

There are a great number of models whose correlation functions are governed by a Painleve 
transcendent; see [M] for a recent review. As pointed out in [34J, in general it is not expected 
that Painleve transcendents arise in correlation functions of models that are exactly solvable 
but which are not free fermion models, such as the six-vertex and eight-vertex models, and 
the XXZ quantum spin chain. There seem to be exceptions that include certain ferromagnetic 
models. The present work goes partly in this direction, as it relates the Painleve V transcendent 
to the six- vertex model with domain wall boundary conditions in the ferroelectric phase, and 
in particular away from the free fermion line of the model. In contrast, two-dimensional Yang- 
Mills theory does have fermion operator representations |35l ll]: this has some implications for the 
usual matrix model description of the six-vertex model with domain wall boundary conditions 
in the ferroelectric phase [14] . 

Another important and natural appearance of Painleve transcendents is in the reinterpretation 
of two-dimensional quantum gravity in terms of matrix models. This approach led, in the work 
of Brezin and Kazakov, Douglas and Shenker, and Gross and Migdal [36, |371 [38], to exact 
solutions of two-dimensional gravity coupled with matter fields. A relationship between the 
free energy in the double-scaling limit of the multicritical matrix models and solutions of the 
Painleve I equation was a crucial result. The role of the discrete Painleve transcendents and the 
tau-function of an isomonodromic deformation in two-dimensional quantum gravity was later 
studied in further detail in [39] [iOl [il] . 

The double scaling limit of the full heat kernel expansion for Yang- Mills theory on the sphere 
is related to Painleve II 02]. In this case the double scaling free energy Fyuit) is given by 
Fy^i't) = v{t)'^ /A, where v{t) satisfies the Painleve II equation 



2v" -v^ + tv = Q 
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The same property is shared by the Gross- Witten model [33]; this is another combinatorial 
model of two-dimensional quantum Yang- Mills theory obtained as a one-plaquette model based 
on the Wilson action instead of the heat kernel action in (jl.2p . While Painleve I appears as the 
universality class of two-dimensional gravity, the Painleve II equation describes two-dimensional 
supergravity [Hj. In the Gross- Witten model the relationship with Painleve transcendents goes 
beyond the connection between the specific heat and Painleve II in the double scaling limit, 
and it involves the Painleve III and Painleve V equations as well |45t 135] . The appearance of 
one equation or the other depends on the matrix model quantity being considered and whether 
or not a double scaling limit is taken; we summarize these results in Appendix A. This feature 
is particularly important in our situation, because by considering the works |47[ 138] we can 
establish analogous results for chiral Yang-Mills theory on 5^ and also for the six-vertex model 
with domain wall boundary conditions. 

Both the six-vertex model and the chiral two-dimensional gauge theory are described by 
a discrete random matrix model, the Meixner ensemble, and this has implications for both 
systems. (The ferroelectric phase of the six- vertex model with domain wall boundary conditions 
was already studied using the Meixner ensemble in @9j.) In particular, it will allow us to 
derive a unitary matrix model description for chiral Yang-Mills theory on 5*^ and also to give 
a stochastic interpretation of the partition function, related to the one given in [28j. Another 
novel consequence is that its partition function satisfies the Painleve V equation. These results 
also apply to the six-vertex model in the ferroelectric phase, but in that case one must always 
consider a certain limit of the model, as we explain later on. 

1.5. Gauge theory characterization of phases in the six-vertex model. Following [4:9\ 
[50l [51j we study the asymptotic large behaviour of the free energy of the six-vertex model 
with domain wall boundary conditions for all three of its phases: the ferroelectric phase (which is 
related to chiral Yang-Mills theory on 5^), the disordered phase, and the antiferroelectric phase. 
Recall that the usual large A^ expansion of the partition function Zj\f follows the topological 
expansion 



introduced by 't Hooft [52] and studied later on in the context of matrix models [53\ 154] . In 
particular, in [53] a model with a quartic potential was shown to follow the asymptotic expan- 
sion p.6p . and more recently the topological expansion has been rigorously proved for matrix 
models with polynomial potentials in [55]. These matrix models lead to a determined moment 
problem [56], but matrix models with weaker confining potentials behave very differently. In 
particular, Chern-Simons gauge theory is described by matrix models whose weight function 
leads to an undetermined moment problem [57\ [58] , and indeed the related large A^ expansion 
of closed topological string theory contains terms that do not appear in (jl.6p [59j . 

As we show in the following, the disordered and antiferroelectric phases of the six-vertex model 
with domain wall boundary conditions exhibit a behaviour that goes beyond ()1.6p and that is 
typical of simple gauge theories, such as two-dimensional Yang-Mills theory or Chern-Simons 
theory. In the disordered phase, for example, we find a term of the form /t log A^, which is absent 
in ()1.6p but present in Chern-Simons theory on the three-sphere even in the semiclassical 
limit. The antiferroelectric phase, on the other hand, exhibits an oscillatory behaviour due to 
a term involving a theta-function; we shall see that this behaviour is characteristic of a matrix 
model with a multi-cut solution and appears in the large A^ expansion of Chern-Simons gauge 
theory on a lens space. 



oo 



(1.6) 




3=0 
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The SU{N) Chern-Simons gauge theories that we consider in this paper are equivalent [6^ 
l6H [62l [63] at level A; G Z to (/-deformed Yang- Mills theory on S^, with the identification of the 
string coupling = and 

q := e"^" . 

This equivalence was exploited in the context of certain one-dimensional integrable models 
in [17J. In particular, the disordered phase is closely related to the full (non-chiral) topological 
Yang-Mills theory on S = 5"^; in this limit the two-dimensional gauge theory partition func- 
tion computes the symplectic volume of the corresponding moduli space of flat connections [S] . 
This is to be compared with the physically simpler ferroelectric phase, which is related to the 
chiral sector of the non-topological gauge theory that computes the orbifold Euler characters of 
Hurwitz spaces. 

Fig. [1] depicts the standard phase diagram for the six- vertex model with domain wall boundary 
conditions. Tab. [1] illustrates the analogous behaviour of gauge theory for each of the three 

4 b 




Figure 1. Phase diagram for the six- vertex model with domain wall boundary 
conditions, depicting the phase boundaries separating the ferroelectric (F), dis- 
ordered (D), and antiferroelectric (AF) regimes as functions of the Boltzmann 
weights a and b (with c = 1). The dashed circular arc denotes the free fermion 
line. 

phases. These correspondences with Chern-Simons theory or two-dimensional Yang-Mills theory 
are not exhibited by other vertex models. Note that the only case where the deviation from 



Phase 




Gauge theory 


String interpretation 


Matrix model matches 


F 


qN pN^ 


Chiral 


Branched covers 


Yes 


D 




Topological 


Flat connections 


Large N 


AF 




g-deformed 


Topological strings 


Large TV 



Table 1. Gauge theory characteristics in each phase of the six- vertex model 
depicted in fig. [U 



the usual topological expansion of matrix models with polynomial potentials does not appear 
is in the ferroelectric phase, which is exactly the case where there is a more precise match 
in terms of the associated matrix models. It would be very appealing to relate the different 
gauge theory behaviours corresponding to the different phases in terms of the spontaneous 
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symmetry breaking from a disordered state to an ordered state, for example in the breaking of 
the topological symmetry of the gauge theory of the disordered phase via the appearance of a 
Yang-Mills coupling in the chiral gauge theory of the ferroelectric phase. 

1.6. Outline. The remainder of this paper is organized as follows. In f|2]we describe the parti- 
tion function of chiral Yang-Mills theory on from the perspective of random matrix theory, 
identifying the pertinent classical random matrix ensemble as the Meixner ensemble. We in- 
terpret Zyyi (5*^ , SU{N)^ as a particular case of the normalization of the z-measure, which 
unravels the integrability properties of the chiral sector of the gauge theory; for example, it 
shows that the partition function is related to a tau- function of the Painleve V transcendent. 
Using the expressions for the partition function involving Toeplitz determinants and the Schur 
measure, following [28] we show that there is an alternative matrix model description to the 
usual discrete matrix models [641 B2] based on n x n unitary matrix models; as in [28] there 
is also a natural stochastic interpretation of the chiral gauge theory. In ^ we explain basic 
facts about the homogeneous six-vertex model with domain wall boundary conditions for each 
of its phases in the thermodynamic limit. We obtain the explicit mapping between the partition 
function of the ferroelectric phase and that of chiral Yang-Mills theory on S*^ in the large A'^ 
limit, and describe various features of this correspondence. We also suggest that the general 
behaviour of the ferroelectric phase can be captured by generalized two-dimensional Yang-Mills 
theory, and relate this to the characterization of the partition function as a tau-function of an 
integrable Toda lattice hierarchy. We further argue that the disordered and antiferroelectric 
phases are described by Chern-Simons gauge theory on the three-sphere and the lens space 
L(2, 1) = S^/'Ij2 respectively (equivalently particular g-deformations of Yang-Mills theory on 
5^), and describe various properties of the mappings in this case. In Appendix A we summa- 
rize some technical details of the relationship between the Gross- Witten model and Painleve 
equations. 

Acknowledgments. We thank R. Weston for helpful discussions. The work of RJS was sup- 
ported by grant ST/G000514/1 "String Theory Scotland" from the UK Science and Technology 
Facilities Council. The work of MT has been supported by a Lady Davis fellowship at the He- 
brew University of Jerusalem and by the project "Probabilistic approach to finite and infinite 
dimensional dynamical systems" (PTDC/MAT/104173/2008) at the Universidade de Lisboa. 



2.1. Meixner matrix model. We begin by deriving the discrete matrix model representation 
of the chiral partition function ZyyiC^, SU{N)). The irreducible representations A of SU{N) 
correspond to Young diagrams for which the number of non-zero rows n satisfies the constraint 
n < N. The row lengths are denoted by Aj for all i = 1, . . . ,n and they define a partition 
A = (Ai, . . . , An), i.e. Aj > Aj+i > 0. The quadratic Casimir operator in these variables is 



2. Chiral Yang-Mills theory and the Meixner ensemble 



n n 



C2{\) = NY,h + Y,\, (A. + l 



2i) 



1=1 i=l 



while the representation dimensions are given by 



diT^X = - J] (A. 



n {Xi-i + N)\ ^<i<j<N 



i=l 



where |A| = Yli is the number of boxes of A. 
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If one considers only small representations, whose row lengths are all < N, then the Casimir 
eigenvalue is linearized in the large limit. The discrete matrix model then follows from the 
heat kernel expansion (jl.3p by dropping the constraint on the number of rows and is given by [6l] 

TV 

(2.1) z+^(s, 5c/(7v)) = E n n ' 

n^eN i<j 4=1 

where we have rescaled the coupling constant Qs N ^ gs, which is held fixed for — t- oo, i.e. we 
work in the 't Hooft limit of the gauge theory. Throughout we drop irrelevant overall numerical 
constants. 

The full partition function is given by a discrete Gaussian matrix model ^42j. In contrast, the 
matrix model (|2.ip is an orthogonal polynomial ensemble [65j which appears naturally in the 
study of stochastic last passage models [BB]. Its associated discrete orthogonal polynomials are 
the monic Meixner polynomials which are orthogonal with respect to the family of weights 

(2.2) WMeix(n;i^) = f ^ j ^z" , neN 

parametrized by -fC G N and q G (0,1), and which can be expressed through hypergeometric 
functions as Mm(n;i^') = 2-^1 ( K \1 — qy, the orthogonality relation is 
00 J 
V 6JMcix(n; K) M^(n; K) M^,{n- K) = ^, . 

The associated Meixner ensemble is a discrete Coulomb gas model on N with joint eigenvalue 
probability distribution [661 EZ] 

(2.3) PMeix(ni,...,nAr;i^') = -j^j^j— — J| (n^ - re^)^ WMcix (n^; K) . 

i<i<j<N j=i 

Hence the partition function of the Meixner ensemble Z^'^'^(X) with K = 1 and the identifi- 
cation q = e in ()2.2p coincides with the genus /i = chiral gauge theory partition function 
Z+^{S^,SU{N)) in (12T]). 

One virtue of this matrix model formulation is that the Meixner ensemble is more general and 
defines natural extensions of the chiral gauge theory defined by ()2.ip . For > 1, it perturbs 
the linear potential in (|2.ip by the logarithmic potential 

y^(n,) = X;'iog(^ 

s=l 

and in this sense it can be regarded as a discrete Penner matrix model. We will use some 
instances of this extension later on in the context of generalized two-dimensional Yang-Mills 
theory [681169]. 

2.2. z-measure. The z-measure is a two-parameter family of probability distributions on inte- 
ger partitions A (equivalently Young diagrams) that originally appeared in the harmonic analy- 
sis of the infinite symmetric group [70]. It has been understood in further detail more recently 
in [711 [Z2 • Consider the quantity 

(2.4) Q^'^'^ = (l-g)^^'z?„(a) 

defined for z, z' G C, n G N and g G (0,1) by an n x n Toeplitz determinant 

(2.5) Dn{a) = det [a,_j\ , 

l<i,j<n 
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where am = CTmiz, z' ,q), m & Z are coefficients in the Fourier series expansion of the associated 
symbol function 

oo 
m=—oo 

A representation theory definition of Qn ^ ''^ states that if = z then Qn is the distribution 
function of the first row of the random Young diagram distributed according to a z-measure |7H 
172] , There are a great number of interesting properties associated to this measure; in particular, 
it can be expressed as a Fredholm determinant with a hyper geometric kernel that allows one 
to establish a very concrete and detailed connection with the discrete Painleve V equation. We 
show that a particular case of the Toeplitz determinant (j2.4p describes SU{N) chiral Yang-Mills 
theory on S'^. 

For this, we use Gessel's formula which expresses a certain series in Schur functions in terms 
of a Toeplitz determinant. In terms of the Schur measure parametrized by sets of variables 
X = {xi)i>i and y = (yj)i>i, the unnormalized probability distribution for the number of boxes 
in the first row of a random Young diagram is given by |28] 

(2.6) Vn{x, y) := Yl (2/) = ^^(^) ' 

A: Ai<Ar 

where Sa(x) := detjj {x'^^^"' '') / detij (a^"""') are the Schur polynomials and 

oo 

(2.7) Am = Am{x,y) = Y ci+m{x)ci{y) 

1=0 

with ti{x) the l-th elementary symmetric function. The symbol of the Toeplitz determinant 
DNiA) in daSl) is [28] 

oo 

(2.8) ^(0= Yl ^m(x,?/)r = n + ■ 

r?i=— oo i>l 

This symbol coincides with that of (|2.4p for z = z' = N and Xi = yi = ^Jq for i = 1, . . . ,N, 
giving an expansion of (|2.4p in terms of Schur polynomials as 

(2.9) D^ia)= Y 5x{VQ^---,Vq)'= E (1, . . . , 1)' = Yl (dimA)^ . 

A:Ai<n ^ ^ ^ A:Ai<n A:Ai<n 

We have to take into account that the Schur polynomial is non-zero as long as its number 
of variables (N in this case) is at least equal to the length of the partition. Hence, mak- 
ing the identification q = e again, we find that the chiral gauge theory partition function 
■^YM ('^^ ' SU{N)) from ()2.ip is a particular case of the z-measure 

(2.10) (52 , SUiN)) = D^ia) = (1 - g)"^' Qj^'^'" . 

This connection with the z-measure also follows directly from the identification of the Meixner 
ensemble as the discrete matrix model underlying the chiral sector of Yang- Mills theory on S'^ , 
since choosing z = N and z' = N + K — 1 in the z-measure leads to the Meixner distribu- 
tion ()2.3p |73] . In the ferroelectric phase of the six- vertex model with domain wall boundary 
conditions, this was already pointed out in [35] . 



10 



RICHARD J. SZABO AND MIGUEL TIERZ 



As an application of this result, let us explicitly evaluate -^ym {^"^ ' SU{N)y For this, we use 
the Cauchy identity for the Schur measure [28] 



N 



Then with the specialization xi = ■ ■ ■ = xn = Hi = ■ ■ ■ = Vn = y/Q = e~^°/^ leading to ()2.9p . 
we find 

Z^^{S\SU{N)) = {l-e-^^y''\ 

This is precisely the formula which was computed in |64j in a much more cumbersome way, by 
attempting to derive the system of discrete orthogonal polynomials associated to the matrix 
model; these polynomials are of course the Meixner polynomials. Below we consider various 
other applications of the representation (|2.10p . 

2.3. Double scaling limit. The z-measure contains the Plancherel measure which is the prob- 
ability distribution on Young diagrams that appears in the description of four-dimensional su- 
persymmetric gauge theory in terms of random partitions [74^ [75] ; in this sense the z- measure 
can be regarded as a deformation of the Plancherel measure. Moreover, a certain double scaling 
limit of the z-measure leads to the Poissonized version of the Plancherel measure. This measure 
is in turn directly related to the Gross- Witten model. This identifies a double scaling limit of 
the chiral gauge theory in which it coincides with the one-plaquette reduction of lattice gauge 
theory defined in ()1.2p . 

The scaling can be described as follows [76j. Generalizing (|2.9p we expand (j2.4p into Young 
diagrams as 

(2.11) Q;^'-^'^= , 

A : Ai<n 

which defines the mixed z-measures 



(2.12) M,,,,,(A) = (l-g)-'(^)%l^l n {z + j-z){z'+j 

(«,i)eA 



|A| 



on partitions (with the product over the boxes of A) . The Poissonized version of the Plancherel 
measure A^pianch(A) = (dimA)^/|A| ! is defined by 

/^|A| . Him A \ 2 

A^Planch.e (A) = A^Planch (A) ^ = ^ ( ^ ) , ^ > . 

Then one has 

(2.13) lim M,,,,,g (A) = A^Pianch,0 (A) . 

2 z' q=e 

The Poissonized Plancherel measure is intimately related to the Gross- Witten model [77], 
which is defined by the unitary one-matrix model with partition function 

(2.14) ZGwiO)= [ dU exp(Ve Tr(U + U-^)) . 

Ju{N) ^ ' 

Expanding the integrand in U{N) characters and using standard orthogonality relations gives 
an expansion of (|2.14p into Young diagrams as [U] 

(2.15) ZGwie) = e' Yl ^ 

Planch, 6 

(A) . 

A: Ai<Ar 
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We can give a natural physical meaning to this result by interpreting the scaling (|2.13|) at the 
level of chiral Yang-Mills theory on S^. It is a large N limit, but the limits g — >• and z z' q = 9 
give q = and hence 

(2.16) N^oo, gs = 2log{N/V0 ) . 

This limit is consistent with the fact that the heat kernel representation (jl.3p is a strong coupling 
expansion of the two-dimensional gauge theory. Hence assuming that the scaling (|2.13p can 
be interchanged with the partition sum in (|2.1ip . we find that the chiral partition function 
Zyjyi [S"^ , SU (N)) from (j2.10p in the double scaling limit (|2.16p is equal to the partition function 
()2.14p of the Gross- Witten model in the large limit 

(2.17) hm {S^ , SU{N)) = e Zgw{0) • 

g,=2\og(N/Ve) 

This result implies that the two gauge theories coincide in the N ^ oo limit. It also shows why 
we need to consider the z-measure, even if only a very particular case, in order to understand the 
connection between Zyy^ (5"^, 5'C/(iV)) and Painleve transcendents for any N. The dependence 
on in that case comes from a very particular specification of the two additional complex 
parameters z and z' in the z-measure. This specification of z and z' to integer values enjoys 
additional symmetries described in the work f78j , where it is shown to imply a connection with 
the counting of branched covers of S"^ . This result fits in nicely with the fact that the chiral sector 
of two-dimensional Yang-Mills theory is related to Hurwitz theory [9l |T0] ; generally the large N 
expansion of the chiral partition function Zyy^ (E, SU{N)) can be expressed as a Gross- Taylor 
string series in branched covering maps of the Riemann surface S. Moreover, consideration 
of the general z-measures ()2.12p can be interpreted as generalized two-dimensional Yang-Mills 
theories |681l69j which are more general than those given by the Meixner ensembles with K > 1. 

2.4. Unitary matrix models and tau-functions of Painleve equations. The expression 
()2.10p for the partition function Zyj^(5'^, SU{N)) in terms of a Toeplitz determinant immediately 
leads to a unitary one-matrix model of chiral Yang-Mills theory on S"^. This follows from the 
Heine-Szego identity j2S] 

[ dU det a{U) = Dn{a) 

JU{n) 

which relates a generic nx n Toeplitz determinant to the integral of its symbol over the unitary 
group U (n) . In the case at hand this leads to 

(2.18) (5^ , SU{N)) = [ dU det (l + f^"^)^ det (1 + ^9 U)^ . 

Ju(N) 

Notice that the rank of the gauge group also appears as a parameter that modifies the 
weight function. The use of the more generic parameters z and z' instead of A^ as powers in 
()2.18p should describe a generalized two-dimensional Yang-Mills theory. Another implication of 
that matrix model description of the partition function is an equivalence between such a matrix 
integral and the one given by the Meixner ensemble. In the double scaling limit considered in 
H2.3\ the matrix model (j2.18p converges to the large A^ limit of the Gross- Witten model (j2.14p . 

An immediate consequence of the expression for the chiral partition function in terms of the 
random matrix average (I238|) is that, by the work 08], the partition function Z^y^(S'^, SU{N)) 
is a tau-function of the Painleve V equation. A systematic comparison with the Gross- Witten 
model is also possible. Forrester and Witte found that the two associated matrix integrals are 
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special instances of tau-functions 

(2.19) [iV] (t; /i) = / dU det (Uf") exp ( ^ Vi Tr (U + U-^)) , 

Ju(N) ^ ' 



and 



\N\ {t; /i, z^) = / dU det(l + U)'' det (l + U'^f exp {t Tr {U)) , 

Ju(N) 

for the Painleve III and Painleve V equations respectively [48], in the sense of the Hamiltonian 
formulation of the Painleve equations (see Appendix A for details). In both cases one can 
identify a Backlund transformation that can be used to establish a Toda chain equation for the 
corresponding tau-function sequence. 
We then have 

^GW(^) = T™ [N] (40; 0) 

and 

hm (^2 ^ ^f^(^)) ^ (Q. ^) _ 

Not only can the partition functions of the two Yang-Mills theories be identified with tau- 
function of Painleve equations, but certain averages in the corresponding matrix models also 
have this property. Averages of determinants in random matrix ensembles often appear in the 
random matrix theory approach to quantum chromodynamics (see |79) for a review) and, in 
particular, the matrix integral (j2.19p appeared in [80] describing a QCD partition function in 
the sector of topological charge /i. Hence the partition function given in [80l eq. (9.25)] is a 
tau-function of the Painleve III' equation, a property already exploited in [79] . 

Notice that we have mostly emphasized the connection with Painleve equations, whereas it 
is well-known that there exists other relationship with integrable hierarchies, mostly of the KP 
type [81j. Tau functions in the Hamiltonian formulation of Painleve equations are related to tau 
functions of the Toda lattice [iH], so both results are related and we expect to explore such a 
relationship elsewhere. 

2.5. Corner growth model. The models considered in this paper are naturally related to the 
corner growth model with geometric weights [UB]; this model is described in terms of the Meixner 
ensemble, and hence is intimately related to the chiral sector of Yang- Mills theory on S"^ studied 
in this section and to the six-vertex model with domain wall boundary conditions studied in ^ 
Let oj for G be independent geometric random variables and define 

G(M,N)= max V u)(i,j) , 

(«j)G7r 

where the maximum is taken over all up/right paths tt in W from (1,1) to {M,N). This model 
can be given several probabilistic interpretations - as a randomly growing Young diagram, as a 
totally asymmetric one-dimensional exclusion process, as a certain directed polymer in a random 
environment at zero temperature, or as a kind of first-passage site percolation model. 
We suppose that w (i, j) are distributed according to the probability measure 

V [w (i. A;) = m] = (1 - g) , m G N . 

In this model, Johansson proved that the probability for G{M, N) being smaller than a certain 
value can be expressed as [661 EZ] 

N / 

P[o,M.^,<,].^ n («.-".)^n "■"! >'" 

^"'^^ rum i=l ^ * ^ 

maxi ni<.t+N —1 
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and hence the normahzation constant of the process Zm,n coincides with the partition function 
of the Meixner ensemble Z^'^'^(A') for K = M — N + \. In the symmetric case M = N, this is 
just the chiral partition function of Yang-Mills theory on S"^, i.e. Z^^n = ^ym('^^' SU{N)). This 
result is consistent with the double scaling limit (|2.17|) : if one takes q = a/N'^, then G{N,N) 
converges in distribution to L{a) as N ^ oo, where L (a) denotes the Poissonized version of the 
random variable describing the longest increasing subsequence in a random permutation, whose 
distribution is given by the Gross- Witten model. 

3. Gauge theory descriptions of the six- vertex model 

3.1. Six- vertex model with domain wall boundary conditions. Let us consider now the 
structure of the partition function ()1.5|) for the six-vertex model. The domain wall boundary 
conditions are only defined for square lattices, and demand that the external horizontal arrows 
are outgoing while the external vertical arrows are incoming. In this model one has six types 
of vertices (ai, 02, 62, ci, C2), where the subscript 2 refers to an opposite configuration to that 
of subscript 1. Gonservation laws reduce the weights to the homogeneous case Ui = a, bi = b 
and Ci = c, and the partition function depends only on the two parameters ^ and | due to the 
identity 09] 

Z,v(ai,a2,6i,62,ci,C2) =c^'z^(f,f,^,^,l,l) . 
The parametrization of the Boltzmann weights associated to the vertices is given by 

(3.1) a = sinh(t — 7) , 6 = sinh(t -|- 7) and c = sinh(27) . 

The partition function of this model was expressed, using the earlier work [21], as a determi- 
nant in [221123]. and this was used in |14j to find a matrix model expression. The determinant is 
a Hankel determinant which has a Hermitian matrix model representation (in the same way that 
a Toeplitz determinant has a unitary matrix model representation [28]). The Izergin-Korepin 
determinant formula is 

( sinh (7 - t) sinh (7 + t) ) 

(3.2) Zn = TN 

' n n^- 



n=0 



where tn is the Hankel determinant [T 



TN(t) = det 

l<i,j<N 



-Mt) 



dt*+J- 



with 



sinh(t -|- 7) sinh(t — 7) 

The partition function tn is a tau-function of the Toda chain hierarchy; in particular, it satisfies 
the Toda equation 

TNt'n - {TnY = TN+1 TN-1 ■ 

To characterize the three phases of the model in the thermodynamic limit N ^ 00 one 
introduces the parameter 

(3.3) ^^5^+5^ 

2 a 

The ferroelectric phase is the region A > 1, the antiferroelectric phase is A < —1, and the 
disordered phase corresponds to — 1 < A < 1. The free fermion curve is given by A = in the 
space of parameters. The large N asymptotics of the partition function in each phase can be 
computed by means of the Riemann-Hilbert method |49[ \^U[ EI] and is summarized as follows: 
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(F) The ferroelectric phase is the region where the two parameters satisfy I7I < t, and with 
7 > for any e > as — t- 00 one has 



Zn = CG^ 1 + Ofe 



with C = 1- e"^^, G = e^"* and F = sinh {t + 7). 
(D) For the disordered phase one has |t| < 7 and [HI [23] 

(3.4) Zn = cN^ e^'f {1 + 0{N~^-^)) , e>0, 
where c > is a constant, while 

1 272 ^ . , /7r(cos(2t)-cos(27)) 

K = -. and t = log :^-7T 

12 37r(7r-27) •' ^\ 47 cos (|^) 

(AF) The remaining antiferroelectric phase is the region |t| < 7 of parameter space and the 
— 7- 00 limit of the partition function is given by [51] 

(3.5) ZN = ce4Nil + ^),q)F''" {1 + 0{N-')) , 

where c > is a constant and F is a function on the two-dimensional parameter space 
given by 

^ _ TT sinh(7 - t) sinh(7 + t) e[{0, q) 
~ 27^1(1 + ' 

while the elliptic nome q of the theta-function 

00 

e^{z, g) = 1 + 2 ^ (-1)" 

n=l 

is given by g = e . 

3.2. Ferroelectric phase and chiral Yang-Mills theory. The match between the six-vertex 
model and chiral Yang-Mills theory on S'^ occurs in the ferroelectric phase with A > 1 (and 
< I7I < t). In this regime the partition function r^v can be written as a discrete matrix 
model [14j 

TV 

(3.6) = 2^^' E n - ^i)' « n ^^^^ (27 li) . 

/jSN \<i<j<N i=l 

Instead of studying this model, in |14j the eigenvalues are rescaled as = k/N and then expo- 
nentially small contributions for large are neglected by using sinh (27 n,iV) = e2|')'l'^«^/2 
0(6"^), leading to the matrix model 

(3.7) T^ = (2iV)^'-^ n (^^-^i)' ^"''^^*"'^'^^^"^ + ^(«"'^) ' 

n,eN/N i<j 

which is a Hermitian Meixner ensemble. 

To compare this matrix model directly to (j2.3p . we want to have a normal scaling of the 
eigenvalues, i.e. S N. This implies that the approximation in |14] then requires 7 to be very 
large, if we do not consider the large N limit. Since I7I < i in the ferroelectric phase, we must 
also have t large such that t — I7I > 0. Hence for 7 — )• 00 the partition function of the six-vertex 
model can be mapped to the partition function of chiral SU{N) Yang-Mills theory on via 
the identification 

(3.8) 5. = 2(t-|7l). 
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To write this partition function as a tau- function of the Painleve V equation, we must further 
work in the weak couphng regime; this is the requirement that t — I7I must be very small. In 
this limit, the Boltzmann weights b and c in (j3.ip are very large while a = t — I7I — )• 0; this 
appears to lead to a reduction to a four-vertex model. 

These results also suggest a connection between the limiting expression ()3.7p for the partition 
function in the ferroelectric regime and free fermion systems. Yang-Mills theory on can be 
regarded as a theory of free fermions on a circle [4J . In the full (non-chiral) partition function 
on the sphere is found by computing the scalar product of the wavefunction for N fermions at 
position X = on the circle at time t = with the wavefunction of N fermions at x = at 
time t = T. The chiral part of the Hilbert space in ()1.4p is the left-moving sector of a 
c = 1 conformal field theory This implies that there is a free fermion Fock space 7^+ whose 
partition function coincides with the partition function of the homogeneous six-vertex model 
with domain wall boundary conditions in the ferroelectric phase A > 1 in the limit of |14j . 
Other relationships between free fermion models and the six-vertex model with domain wall 
boundary conditions are known to hold in the disordered phase [H2j. 

3.3. Ferroelectric phase and generalized Yang-Mills theory. It is natural interpret the 
exact result (j3.6p in the context of generalized two-dimensional Yang-Mills theory 168\ [69] , whose 
heat kernel expansion on a Riemann surface T, is given generically by 

2M''{9s,t) = ^ (dimA)^"^'' exp ( - 5s ^ tpCp{X) 
A p>0 

where 

n ^ 

c,(A)^^(A.-.+if n('- ,._,,+j-, ) 

1=1 jj^l 

is the p-th Casimir operator eigenvalue in the representation A. Higher Casimir operators in 
the heat kernel expansion correspond to higher powers of the field strength F in the gauge 
theory action (jl.ip [69]. Indeed, as pointed in [8], the distinctive properties of invariance under 
area-preserving diffeomorphisms, the absence of propagating degrees of freedom, and exact self- 
similarity [5] are not unique to the two-dimensional gauge theory based on the Yang-Mills 
action (jl.ip . In the correspondence between Yang-Mills theory on S"^ and the six- vertex model, 
we are led to consider an additional potential V{1) := log(2 sinh 27 /). At strong coupling 
t > 7 ^ 1, this potential has an expansion 



g -4n7Z 



n 

n>l 

Since any polynomial in the Young tableaux weights Aj can be written as a linear combination 
of Casimir invariants Cp(A), in this regime the partition function of the six-vertex model can 
be mapped to a modification of chiral Yang-Mills theory by infinitely many higher Casimir 
operators. 

Different truncations of this power series expansion could also be studied approximately using 
the Meixner ensemble with more general (not necessarily integer- valued) parameter K. This 
is particularly useful if we recall the integrability structure of the model, i.e. that ()3.6p is 
the tau-function of the Toda chain hierarchy. Tau-functions of the 1-Toda and 2-Toda lattice 
hierarchies have expansions in terms of Schur polynomials that are useful in illustrating how 
the heat kernel expansion of this generalized two-dimensional Yang-Mills theory compares with 
the ordinary one in (jl.3p (or (j2.ip for the chiral case). In our particular case, the squared 
Vandermonde determinant structure implies that the expansion in terms of products of pairs 
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of Schur polynomials in its diagonal form is the one required (see |28] and references therein). 
More precisely one has 



where the coefficients c^^at are Pliicker coordinates of an infinite-dimensional flag manifold and 
are given as determinants. Note that cx^m = 1 in the case of the usual chiral two-dimensional 
Yang-Mills theory; in this way one can interpret the term Yli sinh(27Zj) in (j3.6p geometrically 
in terms of Pliicker coordinates. 

The expression ()3.6p exactly describes the six-vertex model away from the free fermion curve 
A = 0. Using Meixner polynomials, a comparative asymptotic study of the relationship between 
the exact partition function ()3.6p and the limiting partition function ()3.7p is carried out in [39] . 
This raises the question as to whether or not there is any special connection between the general 
theory and a free fermion system; the same question can be applied to the general Meixner 
ensemble ()2.3p . In the context of generalized two-dimensional Yang- Mills theory, this question 
has been answered affirmatively in [83]. In principle one obtains in this way a description of the 
six- vertex model in the thermodynamic limit from the point of view of Hurwitz theory |69[ I83j. 

3.4. Ferroelectric phase and BF-theory. As a warm-up to the gauge theory descriptions of 
the other phases of the six- vertex model with these boundary conditions, we can alternatively 
characterize the gauge theory of the ferroelectric phase as a BF-theory on 5^ with a non-zero 
theta-angle. This is reminiscent of recent studies of topological order in condensed matter sys- 
tems and statistical mechanics [HI], although the connection between a topological gauge theory 
and a statistical mechanics model (or its one-dimensional quantum mechanics counterpart: the 
XXZ model) is of a different nature. As discussed in ^3.31 the quadratic Casimir invariant C2 (A) 
in ()1.3p can be replaced by any other function of partitions A = (Ai, . . . , A„) since its appearance 
is due to the fact that one is trying to reproduce a continuum gauge theory whose action p.ip is 
quadratic in the field-strength. For example, instead of the Yang-Mills term Tr F'^ in the chiral 
sector we can consider a flux term TrF in the non-chiral theory that describes U{N) topological 
Yang- Mills theory in two-dimensions, and is equivalent to a theta-angle. If we consider the limit 
gs = and a non-zero theta-angle 6^0, then the heat kernel expansion ()1.3p for S = 5^ and 
U (N) gauge group reads 



with Ci(A) = {\i — i + 1) the linear Casimir invariant of the U{N) representation A. This 
leads to the Meixner matrix model (|2.ip with the theta-angle as the coupling constant. 

Two-dimensional Yang- Mills theory has an alternative description in terms of a BF-theory [SJ 
3]. For Qs = and 9^0, this is ordinary BF-theory which is a pure topological gauge theory 
with action 



where i? is a u(A^)-valued zero- form and K is the unit area form on T, = S"^. This theory is 
equivalent to the semiclassical limit {k — )• 00) of Chern-Simons gauge theory on x T, with a 
theta-angle. A quick way of seeing this is to notice that in this limit the heat kernel expansion 
()1.3p is just the sum over all dimensions of U{N) representations and, by Verlinde's formula, 
this is the /c — )• 00 limit of the Chern-Simons partition function [EE\. 




X 



Z7 = J](dimA)2 exp(-eCi(A)) 



A 
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3.5. Disordered phase and Chern-Simons theory on S^. The path integral for SU{N) 
Chern-Simons gauge theory on an oriented compact three-manifold M localizes onto a sum 
over contributions from flat connections. When M — )• S is a Seifert fibration, the localized 
gauge theory is equivalent to g-deformed Yang-Mills theory on the base Riemann surface S after 
analytic continuation of the respective coupling constants. On the three-sphere M = S^, there 
is only the trivial flat gauge connection up to isomorphism, and the partition function reads 

Zcs {S^ , SU{N)) = H 2 sin ^ ^ 



This partition function can be written as a matrix integral [SS] 
(3.9) Zc,(S' . SU(N)) = I / n |i e-^?/^». n (2s.nl3 

with Qs = In this description of Chern-Simons theory in terms of matrix models, the 

computations are carried out with q real and it is possible to make contact with Chern-Simons 
theory by simply identifying gg = |57) . The orthogonal polynomials of this matrix ensemble 
are the Stieltjes-Wigert polynomials |57^ I86j. This Hermitian matrix model is equivalent to the 
discrete matrix model \58\ [17] 



2 

n^eZ i<j 

which defines the partition function of the corresponding g-deformed Yang- Mills theory on 5^, 
with the three-sphere regarded as the Hopf fibration S''^ — )• 5^. There is also an N x N unitary 
matrix model equivalent whose orthogonal polynomials are the Rogers-Szego polynomials [86j. 

The corresponding free energy can be suitably expressed in terms of non-perturbative and 
perturbative contributions 

Fes = log Zcs = Fnp + Fp . 
The splitting of the exact Chern-Simons free energy into perturbative and non-perturbative 
pieces has a physical interpretation in Type IIA superstring theory. The non-perturbative con- 
tribution is the logarithm of the measure factor in the path integral, which is not captured by 
Feynman diagrams, and it gives the exact Chern-Simons partition function in the semiclassical 
limit /c — 7- oo p7| eq. (2.8)]. It has the explicit expression 

i^np - log ^ 

The volume of the gauge group is inversely proportional to the Barnes double gamma-function 



(3.10) G2(iV + l)=n 



Af-l 

n! 

n=0 

which has the asymptotic large expansion [87] 

logG2(A^ + l) = ^A^2 i^g^__L iog^r_^Ar2 + ^ log27r + C'(-l) 

(3.11) +y_^LAr2-2<; 

^i^25(2<7-2) 

where C{z) is the Riemann zeta-function and B2g are Bernoulli numbers. The full partition 
function, including the perturbative contribution, is a g-deformation of the Barnes G-function 
(IXTUIl (and hence oi vol{SU{N))). 
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This leads to the exphcit expansion 

(3.12) F,, = ^ (log,. - I) - 1 logA. + C' (-1) + g 5^ . 

where again here gs N ^ Qs is the 't Hooft couphng which is kept fixed for N ^ oo. The 
free energy (|3.12|) has a very precise meaning in terms of closed topological string theory on 
the resolved conifold geometry [HTj. It is given by a sum over Riemann surfaces in the pure 
Coulomb phase, i.e. with a single hole covering the whole worldsheet, and it computes the Euler 
characters of the moduli spaces of genus g Riemann surfaces. See ^7] for equivalent string 
theory and gauge theory interpretations. 

The free energy is thus of the type (|3.4p which describes the thermodynamics of the six-vertex 
model with domain wall boundary conditions in the disordered phase. Comparing ()3.12p and 
(j3.4p . we see that both free energies coincide and the "unconventional" factor N'^, absent in the 
ferroelectric phase, is given by «; = — which yields 

7 = f (^/5-l) . 

Note that < 7 < ^, as required in the disordered phase. The remaining parameter t is then 
fixed by the value of the 't Hooft coupling constant through ^ ( logo's — f ) = /, giving 

3/2 f { ^'-'^ (^^) ~ (^'y) 



" ' 47cos(g) 

Thus depending on \t\ < 7 and hence on the value of 't Hooft coupling gs, the vertex model 
lives at different points on the phase diagram within the disordered phase. This is an especially 
interesting phase, because together with the antiferroelectric phase the important "artic circle" 
phenomena occurs in this phase - there are macroscopically big frozen and random domains in 
typical configurations, separated in the limit —t- 00 by an "arctic curve" . In the ferroelectric 
phase, we already know that the term A'^"^ is absent and that there is a correspondence with 
two-dimensional Yang-Mills theory. The critical line in the phase diagram of fig. [1] between the 
ferroelectric and disordered phases has k = j, while the critical line between the disordered 
and antiferroelectric phases has k = j2 which is the right factor but with the wrong sign. Note 
that the error term in (|3.4p implies that there is no contribution to the free energy, which is 
consistent with the large expansion (j3.12p . 

The Barnes G-function ()3.10p also appears in the denominator of the Izergin-Korepin formula 
(|3.2p . It further arises as the partition function of the Gaussian matrix model 

■ i=i i<j 

Thus the semi-classical Chern-Simons gauge theory on is essentially a Gaussian matrix model, 
which matches the leading behaviour of the six-vertex model. 

On the other hand, the partition function of the disordered phase is not quite a Gaussian 
matrix model, which would provide the exact match with the non-perturbative Chern-Simons 
partition function, but rather it comes from a fiuctuating model. The partition function in the 
disordered phase is described by a continuous matrix model yu 



1=1 2 j<j 

This ensemble can be studied with the monic Meixner-Pollaczek orthogonal polynomials |89^ l90j: 
it should be compared with the matrix model (|3.9p . The full partition function is given by 
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On the free fermion line that crosses both the antiferroelectric and disordered regions in fig. [H 
characterized by the parameters 7 = f and |t| < j, the full partition function is just [90] 



1 7=7r/4 



This means that the matrix model 



is essentially the square of the Gaussian matrix model. Notice that the weight function of this 
matrix model can also be written as oo (A) = e*^/2cosh j vr. 

In [50j the matrix model (|3.13p is also studied with the rescaling \i = N fn/j giving 

-'^ i=i i<j 

where 

^^('^^ = -C^ - ]7 l°g ( sinhAT;.^ 

with ^ = Nt/'j. In the thermodynamic limit, limAr_!.oo VAr(/u) = = — + this is a 

linear confining potential, which generally in a one-matrix model lies on the border between 
a determined and an undetermined moment problem. This yields a very different behaviour 
compared to standard matrix models; nevertheless, it is still a determined moment and hence at 
least in principle it could be compared to the Gaussian matrix model. See [89j for more details 
and a way to compute in such an ensemble with orthogonal polynomials. Finally, let us point out 
that on the critical line between the ferroelectric and disordered phases the partition function 
exhibits a novel sub-leading fractional behaviour N^^'^ [91], which is again different from the 
other phases. Although the precise gauge theory interpretation of this critical line is not yet 
clear, it is known that certain non-Gaussian matrix models with potentials similar to that of 
(j3.13p and with a leading N'^/'^ scaling arise from localization of superconformal Chern-Simons 
theories [92]. 



3.6. Antiferroelectric phase and Chern-Simons theory on S'^/Z,2. Let us finally consider 
the thermodynamic partition function (|3.5p in the antiferroelectric region. As mentioned earlier, 
this type of behaviour does not correspond to the standard topological expansion of matrix 
models. It has become familiar more recently in the study of matrix models with multi-cut 
domain solutions [93tl94j. Indeed, the solution in [51] involves a discrete matrix model; this is the 
discrete counterpart of the continuous matrix model that describes the disordered phase, which 
has a two-cut solution. Intuitively, the oscillatory behaviour can be related to the possibility 
of eigenvalue tunneling from one cut to another [93]; in the case of antiferroelectric order, with 
pairs of opposing dipoles on the square lattice, this tunneling could have an analog in terms 
of flipping the polarization of one of the dipoles. The problem of multi-cut solutions of matrix 
models has been studied with more detail in [94], where a more general asymptotic formula is 
obtained including oscillations to all orders (coming from a theta- function and its derivatives). 
As shown in [95] , and consistently with rigorous results in the case of the Hermitian one- matrix 
model, the oscillatory terms can be resummed in terms of a single theta-function. Various 
applications of this kind of behaviour in the context of topological string theory are discussed 
in [91 Eg. 

In the context of our gauge theory descriptions of the phases of the six-vertex model with do- 
main wall boundary conditions in terms of two-dimensional Yang-Mills theory and semiclassical 
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Chern-Simons theory, the most precise characterization appears to be in terms of Chern-Simons 
theory on the lens space L{2,1) = S'^/Z,2. In this case the isomorphism classes of flat SU{N) 
gauge connections are labelled by A^-component vectors {pi, . . . ,pn) with pi = 0, 1. The parti- 
tion function is given by 

\ ' Pi=0,l 1=1 

^ n 2sin— ^((/c + iV)(p,-p,)+j-i) . 

l<i<j<N ^ ' 

It can be rearranged into a sum over contributions from the flat connections 

N-l 



Zcs = Yl ^cs 



j=0 

which can each be expressed as a matrix integral [59 



.2 



(3.14) = rj—^, / n — e- - n ( 2sinh: 

where pi = ^^^^ — i is the Weyl vector of SU (N), and (pi, . . . ,pn) is any vector with j entries 
equal to and N — j entries equal to 1. This partition function is related to that of g-deformed 
Yang-Mills theory on S'^ given as 

(3.15) ZgYM{S\ SU{N)) = e-f^S»"? ( 2 sinh | (n, - n,) ' ' 



I ^ Slim 

n^eZ i<j 



via application of the Poisson resummation formula to express ()3.15p as the instanton expan- 
sion [HTl [631 



N-l 

(3.16) Z,YM = ^3(0, e-^-'/^^ye,{^, ,-s^V9s^^-^ ^ 

i=o 

where 

00 

e3{z,q) = 1 + 2 Y 9"'^^ cos (n z) . 

n=l 

The contribution ()3.14p to the Chern-Simons partition function is that of a two-cut matrix 
model. Its large N expansion is studied in [96j. In the "Chern-Simons phase" wherein gg = 
is imaginary, the dominant saddle-point contribution to the two-cut free energy for N large and 
even comes from the "symmetric" flat connection labelled by j = to leading orders in 

this case it is precisely of the form (|3.5p with ^96^ eq. (4.18)] 

Fcs = J,Fo(f ,f)+log^3(^, e^o') , 

where once again Qs has now been rescaled to the 't Hooft coupling constant, and the explicit 
expansions of the genus and 1 free energies around Qs = ^ are given by 

F,{t,t) = ^+log(-4)t^ + t2(logt-^)+^-^ + 0(t«), 
Fi(t,t) = 2C'(-l) + ^logiV + ^ + 0(t4). 
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The gauge coupling gg is related to the weight parameters of the six- vertex model by 

TP" ^'^ A 

Fq = and gs - 



7 2N{t + -f) + l' 

Since this correspondence is made in the 't Hooft limit, it is again valid in the semiclassical limit 
A; — )• oo of the Chern-Simons gauge theory. Thus while the disordered phase is described by 
semiclassical Chern-Simons theory on 5'^, here it is described by the same gauge theory on the 
three-manifold S^/Z2- 

In the "Yang-Mills phase" where gs is real, the large N expansion for Chern-Simons theory on 
S'^/Z2 is equal to that on [961 sq. (4.11)], due to exponential suppression of contributions from 
non-trivial flat connections. It is tempting to relate this behaviour to the spontaneous symmetry 
breaking from disordered states (with gs S M) to ordered antiferroelectric states (with gg G i M) 
where electric dipoles are alternatingly opposite in each sublattice; the appearance of a local 
electric field (with overall zero polarization) can then be attributed to the contributions from 
non-trivial instanton sectors. Note, however, that the proper gauge theory description in this 
case should be expressed through the g-deformed Yang-Mills theory ()3.15p . whose instanton 
expansion (j3.16|) can still provide the appropriate oscillatory behaviour of the antiferroelectric 
phase. 

The matrix model that describes Chern-Simons theory on the lens space L{2, 1) = S"^/Z,2 is a 
two-cut matrix model, which is equivalent to g-deformed two-dimensional Yang- Mills theory with 
the discrete matrix model description ()3.15p . The matrix model describing the antiferroelectric 
phase of the six-vertex model is also discrete and is given by the Hankel determinant [51] 

l<i<j<N i=l 

The fact that the Vandermonde determinant here is not the g-deformed one of (jS.lSp signals 
that this phase describes the semiclassical ? — ?• 1 limit. The continuous version of the very same 
matrix model characterized the disordered phase in ^3.51 and the discussion there surrounding 
the linear (instead of quadratic) confining potential applies here as well. 

Appendix A. Painleve equations and the Gross- Witten model 

The six Painleve equations are non-linear differential equations whose solutions, the Painleve 
transcendents, can be thought of as nonlinear analogs of the classical special functions \97\ 198]. 
The general solutions of the equations are transcendental since they cannot be expressed in 
terms of known elementary functions and, consequently, a new transcendental function has to 
be introduced in order to describe their solutions. The Painleve III equation is 

, . , d^w 1 / do; \ 2 1 dw auj"^ + B . 6 

whereas the Painleve III' equation that appears in e.g. [H] is 

(A 2) (Py ^ dy Y I dy ^ ay'^ ^ 13 ^ jy^ ^ 6 

dx^ y \ dx J X dx 2x'^ 2x x'^ y 

Eq. ()A.2p follows from eq. (|A.ip by substituting oj (z) = y{x)/ y/x with x = \ z^. 

The connection between the Painleve equations and the random matrix averages arising 
in two-dimensional Yang-Mills theory and the six-vertex model does not directly involve the 
equations (jA.ip or (|A.2p . but rather involves the Hamiltonian formulation of the Painleve equa- 
tions |971 [99] . In this formulation, corresponding to every Painleve equation there is a Hamilton- 
ian H (p, q) such that, by eliminating the momentum p in the Hamilton equations, the Painleve 
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equation in q follows. One actually has a family of Hamiltonians Hn, indexed by a parameter n, 
which in the correspondence with a random matrix average is identified as the dimension of the 
matrices. The tau-function is then introduced as a function of an independent variable t and 
the parameters by 

Hn = log Tn ■ 

at 

Thus the Painleve equations themselves are not encountered directly, but rather their cj-form 
due to Jimbo, Miwa and Okamoto. In particular, the cr-form of Painleve III' reads 



(^<ii')^ - ^1 ^2 «r)^ + cTm' {^w - 1) (^iir " *<i') 



f 1 - V2 



64 ' 

while the u-form of Painleve V is 

2 ^ 

(A.3) (t 4)^ - (av - t ctV + 2 {a^f + (i/q + ui + 1^2 + 1^3) (Ty) + 4 JJ {''i + f^V) = ■ 

The Painleve equations in u-form follow directly from the Hamiltonian formalism; the Hamilton- 
ian itself satisfies a differential equation and the cj-form of the Painleve equation is the equation 
satisfied by an auxilliary Hamiltonian, defined in terms of the original Hamiltonian and the 
parameters. For example, in the case of the third Painleve equation one has 

The Gross- Witten model is related to the Painleve II equation, an important result in the 
physics literature |43j that is now known to be a particular case of a more general mathematical 
relationship between the matrix model and the theory of Painleve transcendents |45l H6] . We 
have seen in ^2.41 that the partition function of the Gross- Witten model is a tau-function of 
the Painleve HI' equation. In the mathematics literature, the Gross- Witten model ()2.14p is 
studied as a Toeplitz determinant D]\[{t) = deti<ij<N [fi-j] with symbol f (z) = e*^^^^ \ 
where t = y/O. The usual expression for the partition function in terms of the normalization 
coefficients /i„ of the associated orthogonal polynomials is 



N-l 

1, 



(A.4) DN{t) =11 K 

n=0 

If we set = ^ ~ = hn+i/hn, then (t) satisfies a variant of the Painleve V equation 
The relationship with the partition function is given by 



DN{t) = ex-p (^A J ds log (t/s) s^N (s) 



As shown in J45j, eq. (|A.5p is a Painleve V equation that, using the work of Okamoto [99j, is 
known to be reducible to Painleve HI. In |46) it is also shown that if one introduces the quantities 
Wn = Un/Un-1, then Wn Satisfies another special case of the Painleve III equation. 

This connection to Painleve III is related to the well-known result of Periwal and Shevitz [l3] 
by the consideration of the double scaling limit N ^ 00 and t — )■ 00 with t/N finite, which 
reduces Painleve HI to Painleve II [l5] . Recall that Painleve II actually appears as the continuum 
limit of the discrete Painleve II equation satisfied by the coefficients Un as 

{n+l)Un = t {Un-l + Un+l) (l " U^) , 
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a property that readily follows from ()A.4p together with generic recursion properties of the 
normalization coefficients This property is also used in the derivation of the Painleve III 
equation mentioned above [i5] . 

We have seen in this paper that the chiral sector of Yang- Mills theory on S = S'^, based 
on the heat kernel action, is intimately related to the Gross- Witten model which is also a 
two-dimensional lattice gauge theory, but on S = and based on the Wilson action. This 
relationship only holds for N ^ oo. For finite it is possible to obtain a relationship between 
this chiral theory and Painleve transcendents by using mathematical results concerning the 
Toeplitz determinant (j2.5p . instead of that above associated to the Gross- Witten model. Both 
Toeplitz determinants are studied in detail in |47j . 
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